From microstructural features to effective toughness in disordered brittle solids 
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The relevant parameters at the microstructure scale that govern the macroscopic resistance of 
heterogeneous materials are investigated theoretically. Focusing on brittle failure of disordered solids 
modeled as the propagation of a crack front in a random medium, our study reveals two regimes: 
in the collective pinning regime, the critical failure load can be expressed as a function of a few 
parameters only that can be extracted from the material microstructural features; in the individual 
pinning regime, the passage from micro to macroscale is more subtle and the full distribution of local 
toughness is required to be predictive. Such a rationalization of the disorder induced toughening in 
brittle failure open new perspectives for the design of stronger materials. 



Bridging microscale properties of material with their 
effective behavior at the macroscale is a major challenge 
in both pure and applied science. In systems like brit- 
tle solids [Tj, ferromagnets [2j [3j, superconductors [4], 
martensitic solids [5^, etc. impurities or defects present 
at the microstructure scale can produce dramatic effects 
at the macroscopic scale, and sometimes have interest- 
ing benefits; for example large precipitate particles trap 
dislocations in metallic alloys, increasing notably their 
overall strength [6]. 

In this study, we address the challenge of determin- 
ing the effective toughness of brittle solids in presence of 
disorder, like randomly distributed defects or impurities. 
Roux et al. have proposed a homogenization procedure 
that predicts the shape of the distribution of macroscopic 
toughness [3 [8]. However, it does not give any hint on 
the actual value of the effective toughness and its rela- 
tionship with microstructural parameters. More recently, 
this problem has been addressed in the context of thin 
film peeling [9] . It was shown that periodic arrangements 
of strong impurities could dramatically affect the peeling 
front behavior and increase the overall adhesive perfor- 
mance. However, the case of a random distribution of 
heterogeneities remains largely unexplored. 

Our study reveals two regimes, depending on the dis- 
order strength. Strongly disordered materials are in the 
individual pinning regime, and their behavior depends on 
many microscopic parameters. On the other hand, the 
collective pinning regime occurs for weaker disorder; the 
effective toughness then depends on a few measurable 
microscopic parameters and can be predicted. 

A microscopic description of a crack is in general a very 
complex problem involving dissipative mechanisms such 
as damage and non-linear processes resulting from the 
high level of stress near the crack tip. A major simplifi- 
cation occurs for cracks propagating in so-called brittle 
solids where the typical length scale associated with these 
dissipative mechanisms is much smaller than the typical 
size of the microstructural features in the material. Un- 
der this condition, a crack can be described as a front line 
propagating through quenched heterogeneities p!QHT2] , as 
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FIG. 1: (Color online) Representation of a crack front mov- 
ing upwards in a disordered medium: dark red stands for 
high resistance areas, and light blue for low resistance areas. 
The defects have characteristic sizes and in the r and u 
directions, respectively. 

shown in Fig. [l] The crack dynamics is then governed 
by two ingredients: the material elasticity that tends to 
keep the front flat and the impurities that attract or re- 
pel the front. More quantitatively, a crack front ii(r, t) in 
a material with fluctuating local toughness f{r^u) sub- 
mitted to the external loading /ext follows the evolution 
equation 

du , . ^ of u(r't)—u(rA) , , / 
_(r,t) = + -j -^^-l—±^dr' - f{rMr,t)). 

(1) 

The elastic constant c of the line equals the driving force 
/ext- The long range elasticity follows from [13] and also 
describes the elastic energy of a wetting front ^] . Large 
front perturbations would require non linear corrections 
[T5l[T6j, the effect of which will be neglected here. 

For a random toughness field /(r, li), brittle failure be- 
longs to the wide class of disordered elastic systems, that 
also embraces domain walls in magnetic thin films [17] 
or contact lines during wetting of solids by liquids [TSl - 
|2Q] . In these systems, a complex energy landscape with 
many metastable states emerges from the competition 
between disorder and elasticity. If the applied force /ext 
is small, the crack front line remains pinned by the het- 
erogeneities. When the drive exceeds a threshold value 
/c, the line unpins and acquires a non-zero asymptotic 
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velocity [2T]. 

To describe this behavior, it has been particularly 
fruitful to consider the depinning transition as a regu- 
lar critical phenomenon, with the velocity playing the 
role of an order parameter [22]. The analogy with equi- 
librium critical phenomena suggests that close to /c, the 
front displays universal behavior with critical exponents 
and scaling laws that have been extensively investigated 
[23H26] . If the test of these predictions on experimen- 
tal examples has been rather successful [27H29]. the most 
relevant quantity from an applied science perspective is 
the value of depinning threshold /c. Analogously to the 
critical temperature in equilibrium phase transition, the 
value of fc is not universal and depends, to some extent, 
on the microscopic details of the system. The key point 
is how many and what features at the microscopic scale 
contribute in determining the value of the critical force. 

In engineered composites, the mechanical properties of 
each phase and their concentration allow for the quanti- 
tative characterization of the microstructural disorder. 
For natural materials, high resolution microscopy tech- 
niques allow for the calculation of the toughness corre- 
lator C{5r,5u) = {f{r,u)f{r + (5r,^x + 5u)) - {ff that 
will be shown to provide the relevant microscopic pa- 
rameters; (•) denotes the average over space coordinates 
(r, ix). From a theoretical perspective, the material mi- 
crostructure is modeled by rectangular domains of con- 
stant toughness. Their length is in the r direction; 
in the u direction, it is drawn in an exponential distri- 
bution of average On each domain, the disorder is 
/ = (/) + crF, where F is drawn from a symmetric prob- 
ability distribution P{F) of unit variance and zero mean 
value. Since there is no correlation between the disor- 
der of adjacent domains, the toughness correlator follows 
C{5r^5u) = max ^1 — |-,0^ representative 
of the behavior of most materials. The disorder is thus 
completely described by the three parameters a, and 
and the probability distribution P{F). The average 
toughness (/) can be absorbed in the external force and 
does not play any role in the study of Eq. ([T]). 

Using the symmetries of ([T]), or alternatively a dimen- 
sional analysis, we can show that for a given disorder 
distribution P{F)^ the dimensionless disorder parameter 



allows us to write explicitly the observables dependence 
on the microscopic parameters: 

/, = af;(S,P(F)). (4) 

To perform numerical simulations of Eq. ([T]), the sys- 
tem is discretized in the r direction with a step and 
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FIG. 2: (Color online) Stable front configurations for bivalued 
and rare disorder distributions and two values of the disorder 
parameter S. The fronts propagate upwards and the u coor- 
dinate has been shifted. 

put on a strip of finite width L, with periodic boundary 
conditions. The external force is replaced by a parabolic 
driving centered at w and of curvature hz. We compute 
the first stable configuration Ust{r; w) reached by a front 
initially prepared at u{r) = — oo for 1000 samples. This 
can be done efficiently using the algorithm proposed in 
j3Q] , based on the Middleton theorems [31] . 

In each stable configuration, we compute the roughness 
of the line, defined by 

B{5r) = {[u{r) - u{r ^ 8r)]^) , (5) 

where ~ is the average over disorder, and the pinning 
force /st(^) = ^['^ — {u^ti^f^'^w))] where (•) is the average 
over the position r. The critical force is then obtained 
with the formula /c = /st(^) + ci ^^^^^^^^^ + ' ' ' ^ where 
^/st is the standard deviation of the pinning force distri- 
bution, used in the limit ^ and keeping 1. In 
this limit, \fnL5fJ^ - n^-^ [32 . 

We used four different disorder distributions P{F): bi- 
valued, "rare" (with density n = 0.1) [34], Gaussian, and 
exponential, is varied from 0.001 to 30 and a from 1 
to 8, and we set c = 1 and = 1. The parabolic driv- 
ing force is used as a convenient intermediary tool for 
the simulations, and has no noticeable effect in the limit 

Examples of stable front configurations are shown in 
Fig. [2j for bivalued and rare disorder distributions. At 
strong disorder (S = 1), the configurations look similar 
but the front amplitude is larger for the rare distribu- 
tion. On the other hand, at weak disorder (S = 0.1), the 
configurations look very different but, contrary to the 
previous case, their amplitudes are the same. 

The corresponding roughness functions are given in 
Fig. |3j For a strong disorder, the roughness is much 
higher for a rare disorder than for a bivalued disorder: 
this agrees with the direct observation of the fronts. For 
a weak disorder, the roughness does not seem to depend 
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FIG. 3: (Color online) Crack front roughness: comparison be- 
tween a bivalued (open symbols) and a rare distribution (full 
symbols), for weak (circles) and strong (squares) disorder. 
The solid line is the prediction using Larkin's approximation 
and the dashed line the long distance approximation [Eq. ^\\. 



on the disorder distribution P{F). The difference in the 
front shape is not seen in the roughness; this is not sur- 
prising since the roughness, as a two-point correlation 
function, cannot give a complete account of the front 
shape. 

The disorder induced toughening (or net critical force) 
/c = /c — (/) is plotted versus the disorder parameter 
H on Fig. [4j and two regimes can be distinguished. At 
weak disorder (S < 1), we observe a beautiful collapse 
between different disorder distributions: this is the col- 
lective pinning regime. The net critical force follows the 
phenomenological law 



/c 



= (tS; 



(6) 



that will be justified later by simple arguments. On the 
contrary, at strong disorder (S > 1), the critical force 
depends strongly on the disorder distribution: we are in 
the individual pinning regime. 

In the limit S ^ oo, that can be interpreted as the 
limit c ^ 0, the front is softer and has access to more 
and more disorder realizations. It may thus "choose" 
the more pinning one. As a result, when the disorder is 
bounded, the net critical force gets close to the maximum 
of the distribution. It is indeed what is observed for the 
bivalued distribution, whose maximum is 1, and for the 
"rare" distribution, whose maximum is 3.16. If 

the disorder is not bounded, the critical force is likely to 
diverge, as observed in the simulations. The fatter the 
distribution tail is, the faster the divergence is expected, 
that is also what is observed: the divergence is faster for 
the exponential distribution than for the Gaussian one. 

We now show that simple arguments adapted from the 
Larkin and Ovchinnikov study of vortex pinning in su- 
perconductors ^33^ allow for an intuitive interpretation of 



FIG. 4: (Color online) Disorder induced toughening: The net 
failure load /c is shown for different disorder distributions as 
a function of the disorder parameter S. The line shows the 
prediction of Eq. (10). 



our results. The main difficulty in dealing with Eq. ([T]) 
is its non linearity coming from the disorder term. The 
first Larkin assumption is that this difficulty can be cir- 
cumvented at short distances^ where the crack front does 
not see that the disorder correlation length is finite. 
The Larkin model is thus defined as the limit = oo 
(i.e. S = 0), that amounts of removing the u depen- 
dence of the disorder: f{r^u) fir). The propagation 
equation ([T]) is thus linear and admits a unique stable so- 
lution, that can be determined analytically. This gives, 
after an average over disorder, an exact expression for 
the roughness. In the limit L ^ oo and ^ oo, and at 
long distances r it takes the form 



B{Sr) 



c 



6r 



(7) 



The exact expression and this approximation do not de- 
pend on the disorder distribution but only on its second 
moment <j^. 

The exact Larkin roughness and its limiting law ^ are 
compared to the simulations on Fig. [3) The agreement 
between the weak disorder simulations and the Larkin 
prediction is strikingly good. The long distance approxi- 
mation is also rather good. On the contrary, strong disor- 
der results are far from the prediction; that was expected 
since the prediction corresponds to H = 0. 

We can address the question of the validity of the 
Larkin regime: what does the assumption of "short dis- 
tances" mean? The answer involves the roughness that 
depicts the amplitude of the front perturbations at dif- 
ferent scales. If these perturbations are smaller than the 
correlation length the Larkin model should represents 
correctly the behaviour of the line; on the contrary, when 
they are larger than a different behaviour is expected. 
This defines the length in the r direction up to which the 
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Larkin approximation is relevant, the so-called Larkin 
length Lc, by 



(8) 



If the Larkin length is larger than the correlation length 
^r, its expression can be derived from Eq. 



air 



S2' 



(9) 



The domains of size Lc that behave as if was infinite 
are called Larkin domains. 

To predict the external force required to unpin the 
line, we compute the average resistance seen by a do- 
main of size Lr in the Larkin model. This resistance 
has a mean (/) and its standard deviation depends on 
L^: if Lr < ^r, the domain sees only one defect and the 
standard deviation is a; if Lr > ^r, the front averages 
the local toughness over Lr/£,r uncorrelated defects and 
the standard deviation is cry'^^/Lr' This is valid up to 
Lr = Lc^ where the Larkin model breaks down. Beyond 
that scale, the important wandering of the crack front 
results in correlations of the local toughness actually ex- 
perienced by the system. The second Larkin assumption 
is that the critical force is given by the typical toughness 
seen by a Larkin domain; this leads to two cases: 

• if Lc > (or S < 1), a Larkin domain sees several 
defects: this is the collective pinning regime, where 
the net critical force is /c = G^J^r|L^. 

• if Lc < (or S > 1), a Larkin domain sees only 
one defect: we are in the individual pinning regime. 
The net critical force is now /c = a. 

Using the long distance expression of the Larkin length 
([7|), we can write explicitly the net critical load. 



( Cr^^r V -f V / 1 

= crL II L = —— < 1, 



if S 



(10) 



> L 



These expressions are compared to the simulations in 
Fig. [4] for the four different disorder distributions. In 
the collective pinning regime, the prevision captures the 
results of the simulations: this confirms the validity of 
Eq. ([6|, irrespective of the underlying local toughness 
distribution. The effective toughness, that is a priori a 
non-universal quantity, is shown here to be remarkably 
robust. 



In the individual pinning regime, the prediction (10) 



does not account for all the investigated materials, but 
provides a lower-bound of the effective toughness, as is 
seen in Fig. |4] Indeed, the pinning process filters the 
local toughness distribution towards strong defects of re- 
sistance / > cr. 



Since the elastic constant of the front is given by the 
critical force, Eq. (|6| for the collective regime is actually 
an implicit equation for /c, whose solution is 



/c 




(11) 



It follows that the disorder parameter satisfies S < 
yirjiu] this notably implies that the individual pin- 
ning regime is only accessible to anisotropic materials. 
The expression above represents a powerful tool to de- 
sign the microstructure of composites in order to achieve 
improved failure properties. 

To conclude, we have shown that the interaction be- 
tween the front of a brittle crack and a disordered mi- 
crostructure is governed by one dimensionless parame- 
ter and that two distinct regimes can be identified. In 
the collective pinning regime that occurs at low disorder 
strength, many impurities act together to pin the crack 
front. These impurities are averaged and give rise to an 
effective toughness given by Eq. (|6|, that depends on a 
few parameters measurable from the material microstruc- 
ture. On the other hand, at high disorder strength, the 
front is pinned by strong individual defects: this is the 
individual pinning regime. The basic quantities studied 
here are shown to depend on many parameters, such as 
the disorder distribution. This has a negative effect on 
the predictability of systems falling in this regime, but it 
might be used advantageously for the design of compos- 
ites with improved toughness. Notably, our study shows 
that a few but strong impurities produce a larger tough- 
ening than a large number of weak ones. 
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